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We present a new method for detecting the missing baryons by generating a template for the 
kinematic Sunyaev-Zel'dovich effect. The template is computed from the product of a reconstructed 
velocity field with a galaxy field; we find that the combination of a galaxy redshift survey such as 
SDSS and a CMB survey such as ACT and PLANCK can detect the kSZ, and thus the ionized gas, 
at significant signal-to-noise. Unlike other techniques that look for hot gas or metals, this approach 
directly detects the electrons in the IGM through their signature on the CMB. The estimated 
signal-to-noise for detecting the galaxy-momentum kSZ cross-correlation is 4, 9, and 12 for ACT 
(with survey area of 2000 deg^) with SDSS-DR4, SDSS3 and ADEPT respectively The estimated 
signal-to-noise for PLANCK with SDSS-DR4, SDSS3 and ADEPT is 11, 23, and 32. Our method 
provides a new mean for determining properties of the ionized gas in the Universe. We provide 
galaxy momentum templates constructed from Sloan Digital Sky Survey online at our website at 
http://www.astro.princeton.edu/~shirley/SZ/SZ.html, The predicted correlation coefficients 
are provided along with the momentum maps. One can download the momentum templates and 
cross-correlate directly with CMB maps from ACT and PLANCK to detect the missing baryons. 



Where are the baryons? Astronomers can measure the 
baryon abundance three minutes after the big bang us- 
ing Deuterium [1], 300,000 years after the big bang using 
the angular power spectrum of the Cosmic Microwave 
Background [4], and at redshift three using the Lyman- 
a forest H. However, today, most of the baryons are 
missing [J|. Too hot to show up in Lyman-a absorp- 
tion studies, too cool to find in spectral distortions of the 
cosmic microwave background, and hidden in large-scale 
structures of the cosmic web of too low a density to ap- 
pear in the X-ray, finding this missing matter is one of 
the open challenges in cosmology. 

The baryons within clusters - shock heated during viri- 
alization of the structure - have been the subject of ex- 
tensive study with X-ray telescopes [1, J, 7]. However, 
as X-ray emission strength is proportional to the square 
of electron density, these studies can only detect baryons 
in the very center of clusters. There have been numerous 
attempts to make X-ray observations of the filamentary 
gas Q, but filamentary gas seems to elude X-ray detec- 
tion. 
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Baryons within clusters are now being detected at sig- 
nificant confidence levels through their distortion of the 
primordial black-body spectrum, both in observations 
targeted towards known clusters |9| , in cross-correlation 
with the WMAP data 0, [O, [3i and now in blind 
surveys as the latest detectors come online [ll]. This 
effect, known as thermal Sunyaev-Zel'dovich (hereafter 
tSZ, is, at the arcminute scales to be probed by 

current and upcoming experiments such as the Atacama 
Cosmology Telescope (ACT), the South Pole Telescope 
(SPT) and the Planck satellite, the strongest signal (in 
two-point correlation) on the microwave sky. Measure- 
ments of the tSZ effect can detect baryons in clusters. 
The tSZ, however, depends for its strength on both the 
density and the temperature of the gas, and so an exper- 
iment that may detect the gas at the center of a cluster, 
where temperatures can reach over 10'' K, will be unable 
to find at all the gas associated with the far less dense, 
and far colder, material that lies in the filaments of the 
cosmic web . 

Even when a survey reaches sensitivities of 2 
fiK arcmin"^, identifiable tSZ sources are all objects with 
mass greater than 2 x [l^, and most of the con- 

tribution to the power spectrum comes from virialized 
gas in groups just below the detection mass p7j . 
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The kinetic Sunyaev-Zel'dovich effect [3, hereafter 
kSZ], in contrast to the tSZ, depends not on the ther- 
mal motions of the gas, but rather on the bulk velocity 
of the structures the gas inhabits. In particular, the flows 
of matter toward overdensities Doppler-shift CMB pho- 
tons to produce anisotropics. The fractional change in 
temperature due to kSZ, 8 = AT/Tcmfc, is 

rno 

e(n) = - / dr] g{f])n-p{hr],7]), (1) 

where n is the unit vector pointing away from the ob- 
server and the momentum field, p, is defined as 

p(n?7, 77) = [1 -I- 5b{nf])]v],{n'q, 77), (2) 

and Sb is the baryon overdensity, and the Vb is the baryon 
velocity. The visibility function g is defined as 

g{rj) = x^thH + zfe--(^-\ (3) 

and Xe is the ionization fraction, th is the Thompson- 
scattering optical depth to the Hubble distance today; r, 
in the reionized epoch, is 

r{z) = - + + zf - 1] (4) 

and, finally, rj is the comoving distance defined in units 
of Hubble distance: 

Baryons in the cosmic web have sufficient velocity and 
column density to produce a detectable kSZ signal p^ . 
Indeed, at arcminute scales, the kSZ is expected to be the 
strongest anisotropy after subtraction of the tSZ using 
frequency information (loj . Here we demonstrate a new 
method for extracting the kSZ that is not subject to the 
systematics one might expect when trying to see it in the 
two-point CMB auto-correlation. 

In particular, we consider the idea of cross-correlating 
the CMB sky with a kSZ template constructed by pro- 
jecting down the line-of-sight momentum field of large 
scale structure. This momentum field is reconstructed 
from the inferred three-dimensional dark matter distri- 
bution. We construct such a template from the latest 
large scale structure data, and make it available to the 
scientific community. 

To generate the momentum field we require knowledge 
of the velocity field. We can rely upon the linearity of 
the large-scale velocity field, v(k), which is related to the 
density field via the continuity equation, 

d\nGaHSm{k)k 
"^'^)=^ rflna P ' 

where a is the scale factor, G is the growth factor at 
late times (proportional to a in the matter dominated 
regime) , 6^ is the matter density fluctuation field and k is 



the comoving wave number. This is strictly valid only on 
larger scales where material remains in the linear regime; 
we shall determine how well it holds in approximation 
and on different scales in the following section. 

The temperature increment for a particular velocity 
estimate depends on both the cosmological parameters 
through Eq. [51 and on the gas column density and thus, 
in conjunction with other observations that can constrain 
cosmological parameters, can unambiguously detect the 
missing baryons in the large, low-density filaments that 
have undetectable tSZ. 

There are two major advantages to this method over 
others that have been proposed. The first is that we 
expect the sign of the line-of-sight velocity to be uncor- 
related with many of the standard systematics, such as 
insufficiently subtracted tSZ, galactic foregrounds, and 
detector and telescope noise. 

Secondly, our method makes near-maximal use of the 
information available in a potential cross-correlation. 
Previous studies of constraining dark energy parame- 
ters by using kSZ information in cross-correlation with 
a CMB signal have "thrown out" some of the signal, 
dropping the phase information that determines the sign 
of the kSZ signal, and information about where velocity 
flows are most likely to be found. 

The idea of using peculiar velocity flows in a galaxy 
survey to do cosmology has been studied in past decades, 
most notably by the use of redshift and luminosity dis- 
tance indicators for very nearby (z < 0.1) objects {2^. 
However, the method we propose here has only become 
technically feasible in recent years with the combination 
of large scale redshift surveys and the current and up- 
coming arcminute scale CMB experiment, and the issues 
and questions we must confront are very different. 

Our paper has four parts: the reconstruction of the 
velocity field (Sec. [l|, the determination of the cross- 
correlation coefficients (Sec.|lT]), signal-to-noise estimates 
for current and forthcoming surveys (Sec. IIII[) and the 
template construction (Sec. IIV[) using actual large scale 
structure data from SDSS. 



I. VELOCITY FIELD RECONSTRUCTION 

We first detail our method of velocity reconstruction, 
which provides us with the velocity field that is accu- 
rate in large range of scales which are relevant to the 
structures that harbor most of the gas in the Universe. 
We also validate our reconstruction method via a N-body 
simulation. 



A. Theory 

The Universe provides us with the redshift space over- 
density field, and there are man y w ays we can construct 
a real space density (see, e.g., [22|)- We take here the 
simplest approach, since the scales we are interested in 
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are not strongly affected by non-linear effects - fingers of 
God, for example, occur on cluster scales, whereas the ve- 
locity field has most power on scales larger than 10 Mpc. 
We directly transform the redshift density field into a real 
space density field using the Friedman equation (with 
WMAP5 parameters), and embed the real space density 
field in a box. We compute the overdensity field in real 
space, convert it to Fourier space, and use the transfor- 
mation described by Eq.[6]to determine the velocity field. 
We then convert this field back into real space, and use 
this representation to study the properties of cosmologi- 
cal velocity flows. 

Our method falls into six steps. 

1. Transform the redshift space overdensity field to 
real space overdensity field. 

2. Fourier transform the overdensity field. 

3. Wiener-filter (as discussed below) the map given 
our knowledge of noise from Poisson statistics and 
non-linearities. 

4. Apply Eq. [6]to compute v(fc). 

5. Inverse Fourier transform v(A;) to find the real space 
velocities. 

6. Transform the real space velocities back to redshift 
space, and find the momentum field by taking the 
product with a density survey (which may or may 
not be the same as the survey used to determine 
the velocities.) 

While this method sounds simple, its error properties 
can be hard to estimate. In particular, even if observa- 
tional noise in the density field is uncorrelated (e.g., in 
the case of Poisson noise alone), the noise in the velocity 
field will be correlated from point to point. Formally, 
this can be seen by noting that the Fourier transform of 
1/k is proportional to l/r^, and so in real space com- 
puting the velocity field corresponds to convolving with 
the gravitational force kernel. The induced statistical er- 
rors remain Gaussian - but their spatial correlation is no 
longer a delta function. 

It is also important to recognize that points near the 
edges of the survey will have errors from our lack of 
knowledge of what lies beyond the survey volume. Intu- 
itively, one can imagine the possibility of a "great attrac- 
tor" lying just outside the boundary [l^, [23|, changing 
the velocity field radically from the estimate produced 
by using only the observed material. Developing math- 
ematical techniques to determine the influence of these 
errors on the momentum field is difhcult. We side-step 
the problem in this section by using an N-body simula- 
tion, where we have knowledge of the true velocity field. 



B. Simulation 

We check the validity of our reconstruction method 
with simulations. We first simulate a galaxy field and 



then implement our method of computing the velocity 
field. We then compare our recovered velocity field to 
the true velocity field in the simulation and investigate 
how well the reconstruction has worked. 

Our N-body simulation uses WMAP 5-year parame- 
ters, has volume 256'^ Mpc and contains 256^ par- 
ticles. The linear CDM power spectrum was gener- 
ated using code from the GRAFIC2 [61|| package [23 |. 
GRAFIC2 was then used to generate the initial par- 
ticle conditions, with the modification that the Han- 
ning filter was not used because it suppresses power on 
small scales [1^. Simulations were carried out using the 
TPM [ell code dil with a 256^ mesh and a sphne soft- 
ening length of 20.35 h~^hpc. The initial domain decom- 
position parameters in the TPM code were A = 1.9 and 
B = 8.0 (TPM was modified slightly so that there was no 
lower limit to B when it is reduced at later times, which 
improves the trackin g of low mass halos; for details on 
these parameters see j26j). 

We take this simulation and find its halos using a 
friends-of-friends (FOF) algorithm. After we identify 
the halos, we use the Halo Occupation Distribution 
(HOD) from Bootes field [IJ to populate the halos 
with galaxies. The galaxies were selected in the Bootes 
field by enforcing uniform comoving number density of 
h = 10^^ {mathrmMpc/ h)~^ . We then only take these 
galaxies and compute the overdensity field. This allows 
us to assess the effects of incomplete sampling of the dark 
matter field by the galaxies. 

To quantify the accuracy of the construction, we define 
a velocity reconstruction coefhcient, r(k), as: 



r(fc) 



'(fc)*v''™(fc)) 



^(fc)*V«™"(A;)) ' 



(7) 



where v'''^'^°'^(fc) are the reconstructed velocities, and v*'™ 
those found by the simulation itself. This is a simple way 
for us to gauge how good we can reconstruct velocities 
in the presence of both Poisson noise (from incomplete 
sampling of the dark matter field by galaxies) and break- 
downs of Eq. [nidue to non- linear evolution. 

We apply a Wiener filter to the density field before we 
reconstruct our velocity field, defined as follows: 



W,(fc) = 



(8) 



where P(A:) is the dark matter power spectrum, 6 is a 
(possibly weighted) average bias for the survey, and n is 
the number density of the survey. 

By comparing the top and bottom panels of Fig. [U 
one can see the effect of this filter on the reconstruction; 
while small scale power in the velocity field is filtered 
out, the large scale velocities are traced reasonably well. 
The Weiner filter for generating the density field is sim- 
pler, since there is no additional "reconstruction error" 
as there is for the velocity: 
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FIG. 1: The top panel shows velocity in z-direction for a slice 
in our simulation. The middle panel shows the reconstructed 
velocities of one slice of the simulation when we use all the 
dark matter particles in the box. The bottom panel shows 
the reconstructed velocities of the same slice with Wiener fil- 
tering. The velocities (only in z-direction) are pointing left 
and right of the map in both panel. 



We have two main sources of noise: an inability to 
determine non-linear velocities fields, and Poisson noise 
from galaxies. The effect of the latter can be seen in the 
bottom two panels of Fig. (TJ the middle panel uses all 
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FIG. 2: Smoothed correlation coefficients (l/r(fc)) between 
the reconstructed velocities and the real velocities are shown 
here. At small k, the limited size of the simulation box pro- 
vides us with insufficient information to reconstruct extremely 
large scale velocities - thus the drop off at small k. At large 
k, where non-linear effects kick in, the linear theory we apply 
to reconstruct velocities is not valid anymore. We however 
smoothed out the small k drop off, since this is an artifact of 
the limited simulation box size. We smoothed out the corre- 
lation coefficients at all scale so that we can apply a smooth 
r{k). However, for a large range of scales we have nearly 
perfect reconstruction. 



the particles in the simulation and so is sensitive to non- 
linear breakdowns of Eq. [6l the bottom panel uses only 
'galaxies' to reconstruct the velocities, thus is sensitive 
to both the non-linear breakdowns and also the Poisson 
noise from galaxies. 

As we can see in a plot of r(k) in Fig. [21 the velocity 
field reconstructed from the galaxy density field agrees 
well with both the velocity field reconstructed from the 
matter field and with the full velocity field of the simula- 
tion. There is a significant range of k in the linear regime 
that is well-reconstructed, while the correlation drops off 
at non-linear regime. 

The studies of this section thus suggest that velocity 
reconstruction is a viable method even in the presence 
of the levels of noise we might expect in a contemporary 
large-scale structure survey. We now turn our attention 
to making an analytic estimate of the kSZ signal. 



II. ESTIMATING THE KSZ SIGNAL IN 
^-SPACE 

Having validated our reconstruction method, we now 
turn our attention to predicting the dependence, on cos- 
mological parameters, of the cross-correlation between a 
projected momentum field and the CMB. 
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A. Theory 

The velocity field in the linear regime is a pure gradient 
field: the direction of the velocity is aligned with the 
direction of its A:-mode. Combining Eqs. [5] and [3] in the 
Limber approximation, then, we can see that the source 
of the signal will be primarily from a modulation of this 
velocity field by the density field. This "second-order" 
nature of the kSZ was realized early on [2^. 

There is also the possibility of modulating the velocity 
field by the visibility function: the so called "velocity- 
Doppler" effect that can be significant on large scales. 
However this signal is drawn primarily from redshifts 
around reionization where the derivative of the ioniza- 
tion fraction is large [1^; it does not correlate with a 
low-redshift survey. 

To estimate the momentum-kSZ cross-correlation sig- 
nal, we must first make assumptions about how well our 
various observables trace the underlying dark matter and 
baryon fields. We take, as reasonable Ansatze for the re- 
lationship between baryon, dark matter and galaxy over- 
densities: 



and 



Sb(k) = fimnik) 



(10) 



(11) 



where /(fc) describes the difference in clustering power 
between baryons and dark matter - we take the func- 
tional form of [13] - b{k, z) is the scale and redshift de- 
pendent bias of galaxy overdensity field, and Sm{k) is 
the non-linear cold dark matter overdensity field. We 
take Vgai, our reconstructed velocity field discussed in 
the previous section, to be 



Vgal(k) 



.dlnGai?<5gai(k)k 



din a 



fc2 



(12) 



where, note, we do not divide out by the bias; this allows 
us to keep a handle on how bias estimates will affect our 
final result, and we similarly use /(fc)^gai as our estimate 
of baryon overdensities. 

Once we apply the Wiener filters of Eqs. [8] and [9] to 
our estimates of the baryon velocity field, Eq. [T^ and the 
baryon density field, we may Fourier transform into real 
space and project down to produce the reconstructed, 
projected momentum field, which we call Q. 

We can then ask: how well does Q correlate with Ql 
We here follow and elaborate on the analysis of Ref. [3l| . 
We know that Vgai and Vb are imperfect correlates, and 
assume that r{k), defined in Eq. [71 captures all of these 
effects; this amounts to assuming that baryon velocities 
are governed entirely by gravity and so track the cold 
dark matter - a reasonable assumption on the scales we 
consider. 

We further assume that the effects of Poisson noise, 
non-linear velocity fields, and non-gravitational interac- 
tions are not correlated with the linear field, and that 



differences between dark matter and baryon flows on the 
relevant scales are negligible, so we may write. 



(t.fc(k)t.gal(k')> = {2TTfSD{k-k') 



xb{k, z)r{k) 



dlnGaH 
din a k 



P(fc), 



where P{k) is the non-linear cold dark matter power 
spectrum; note that the correlation here depends also 
on bias. 

We write the angular power spectrum of the cross- 
correlation as 



lQe{l/v), (13) 



where 77 is comoving distance in units of Hubble distance 
as defined in EqlSl G{r]) is the growth factor at r], and 
overdots are with respect to ry. Iqq is the cross-power 
spectrum of the vorticity of the momentum field, which 
we can, after manipulation, write as 



lQeik,z) ^ J dyi J 



47r4 

[bik')fik')P{k')W,,,{k') 
b{k")r{k")P{k")W^{k") - 

4b{k')r{k')P{k')W^,,{k') 

b{k")f{k")P{k")W,{k'% (14) 



k' = kyi, k" = ky2 and 7/2 is equal to \/l — 2/i?;i + y^. 

To compute the signal-to-noise, we will also need the 
auto-correlation of the template, 



{Ih). (15) 



Since we Weiner filter both the density field and the ve- 
locity field in the construction of the template Q, /qq is 
not identical to /qe; we find 

T \ r ^ r ^ (i-M^)fe^ 

Icm{k,z) ^ J dyi J dfi 

[b'{k')P'ik')W^,,ik') 
b\k")P\k")W^^{k") - 

4b'ik')P'{k')Ws,,{k')W,{k') 
2/2 

b\k")P'{k")W,,i{k")W,{k")] (16) 

We have here defined P'{k) as P'(fc) + l/b'^{k)n. 

The "science product" consists both of a template Q, 
defined in this section, and an estimate, (^Q^), analytic^ 
its cross-correlation amplitude with the CMB for a set 
of fiducial choices for the cosmological and galaxy survey 
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parameters. The ratio, R of this estimate to that actually 
measured with an experiment is 



R 



Q0, observed 



^Q0, analytic 
2 



K Ks9cS 



(17) 



where b^s is an "effective" averaged bias, weighted by 
the various kernels 17(77) and N{ri), t/eff is the effective av- 
eraged visibility function within the survey volume that 
has been weighted similarly as the bias. In the case that 
the fiducial parameters chosen are the true ones, R is 
unity. 



Simulations 
Analyticai Calcuiations 




B. Validating Estimates 

Since this is the first time this method has been con- 
structed and applied, we check the validity of our cal- 
culation for through numerical simulations, in the 
following fashion: 

1. We produce the "real" kSZ sky, 8, with an a N- 
body dark matter particles simulation box (as de- 
scribed earlier in Sec. II Bp from z = 0.428 and as- 
suming that the ionized gas traces the dark matter. 

2. We produce the template, Q, with FOF (fricnds- 
of-friends) halos found from the same N-body and 
then use an HOD (from Ref. [27]) to populate the 
halos with their galaxies. 

3. We project these two fields, and take their cross- 



correlation to find C, 



qe.sim 
I 



Note that the highly limited red-shift range of our sim- 
ulations means that we will not reproduce an entire kSZ 
sky - only that produced by baryons in a very limited 
range. 

Fig. [3] compares these two quantities; we find that at 
scales where there are a sufficient number of modes in 
the N-body to do our reconstruction reliably, the analytic 
calculations perform well and without appreciable bias. 

Given the expected signal-to-noise of upcoming mea- 
surements, then, the precision of our analytic formula for 
/qe are more than sufficient for accurate determination 
of gas parameters. 



III. ESTIMATING SIGNAL-TO-NOISE FOR 
THE KSZ CROSS-CORRELATION SIGNAL 

In this section, we evaluate whether current and up- 
coming surveys are capable of using the method devel- 
oped in this paper to detect the kSZ-galaxy momen- 
tum template cross-correlation. Our method requires two 
very different types of observations: a sensitive high res- 
olution CMB map (e.g., PLANCK, ACT or SPT) and 
a large- volume spectroscopic survey (e.g., SDSS, SDSS3, 



FIG. 3: We show the validity of our equations via comparison 
between the expected cross-correlation coefficient (7^®'''"'''^*"^ 
and the calculated cross-correlation coefficient from simula- 
tions CpQ-^'-". 



ADEPT). We estimate the expected signal/noise in the 
combined maps by computing: 



f sky ^{21 



(^CMB 



where /sky is the fraction of sky that is covered by both 
the spectroscopic survey and the high resolution CMB 
survey. Cf^^ is as described in EqUJl while Cf^^ is 
simply the 2-point correlation function of Q (the mo- 
mentum field) with its noise. From the CMB side, we 
have (7p^^ which is the CMB temperature anisotropies 
(which is noise here), and Cf^"^ which is the detector 
noise from the CMB experiment. 

We consider the following two CMB experiments 
specifically: 

1. Atacama Cosmology Telescope (hereafter ACT, 



see 



assume a 
^ with 



for further information): we 
^FWHM of 1.4' and a noise of 26 fiK arcmin 
a survey area of 2000 deg^ based on a straw-man 
proposal. 



2. PLANCK (see [6J| for further information): we 
assume a 6'fwhm of 7.1' and a noise level of 302 
liK arcmin"^ for 75% of the sky, since galactic fore- 
grounds may prove to be hard to subtract from 
parts of the maps. 

For large scale structure surveys, we consider two cur- 
rent surveys and a proposed mission concept: 

1. Sloan Digital Sky Survey (hereafter SDSS): we use 
the SDSS DR4 VAGC LSS sample [13 (which will 
be further described in Sec. IIV[) . The bias of the 
main galaxy sample is found to be 1.2 if we as- 
sume linear bias (see Ref. [2l[). For the purpose 
of signal-to-noise analysis, it is sufficient to assume 
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S/N 


ACT (2000 deg^) 


PLANCK 


SDSS-DR4 


4.1 


11.2 


SDSS3 


8.5 


22.5 


ADEPT 


12.2 


32.2 



TABLE I: The signal-to-noise ratio for the detection of gas- 
momentum-kSZ correlation for different combinations of cur- 
rent and upcoming experiments. 

linear bias, but we will discuss the effects of bias 
later in Sec. El The bias of the spectroscopic LRGs 
is found to be ~ 2 (based on the powerspectrum 
analysis done in Ref. [ssl). We include bias as a 
free parameter throughout the theoretical calcu- 
lation, with one exception. We only employ bias 
from other analysis when we generate filter func- 
tions (VFgai and W^) as defined in Sec. II Bl 

2. Sloan Digital Sky Survey 3 (hereafter SDSS3, see 
[65l | for further information): We assume the avail- 
ability of 10^ spectroscopic Luminous Red Galax- 
ies over a quarter of the sky in the redshift range 
of 0.2 — 0.6. Since SDSS3 plans to spectroscopi- 
cally observe all of the photometric LRGs in SDSS, 
we use the redshift distribution and bias {b ^ 2) 
determined for the photometric LRGs in SDSS as 
described in Ref. [3J|. 

3. Advanced Dark Energy Physics Telescope (here- 
after ADEPT, see [66] for further information): we 
assume the availability of 10^ galaxies over 28600 
deg^ from 1 < z < 2. We assume a bias of 1.5 and 
a uniform distribution in comoving volume from 
1 < z < 2 here. Since these galaxies are Lyman- 
alpha emitters, we take the bias from studies by 
Ref. [3^ which suggests a bias of ^ 1.5. 

For the S/N calculation, we also assume an ionization 
fraction of 1 {xe = 1) and full hydrogen ionization. 

IV. TEMPLATES 

Finally, we illustrate our method by constructing a 
galaxy momentum template form the Sloan Digital Sky 
Survey. This template is a prediction for the AGT and 
PLANCK experimental specifications and this section 
demonstrates how to develop templates for other surveys. 

A. Data from Sloan Digital Sky Survey 

The Sloan Digital Sky Survey has acquired uqriz CCD 
images of 10" deg2 of the high-latitude sky [301. A ded- 
icated 2.5m telescope [13, [3 Apache Point Observa- 
tory images the sky in photometric conditions [3^ in five 
bands [io, HH using a drift-scanning, mosaic CCD cam- 
era [131 • AH the data processing are done by completely 



automated pipelines, including astrometry, source iden- 
tification, photometry (4^ . [43|, calibration [11,14^, spec- 
troscopic target selection [i^, > ^.nd spectroscopic 
fiber placement [i^ . The SDSS is well underway, and has 
produced seven major releases fsol. IsTI. Is^. Issl [54. Issl. l56t . 

In this paper, we utilize mainly the SDSS DR4 VAGC 
(^32^1) sample and also the SDSS spectroscopic Luminous 
Red Galaxies (hereafter LRG, [1^), since these two sam- 
ples are the largest spectroscopic samples publicly avail- 
able and have near-uniform completeness over a large 
area of sky. The spectroscopic LRG sample [46| , includes 
area beyond Data Release 4 (DR4) . The total area cover- 
age for this spectroscopic sample is 5154 square degrees, 
as available in the NYU Value Added Galaxy Catalog 
(VAGC [3^) at the time of the preparation of our project. 

B. Constructing Templates 

We construct the templates by the following steps: 

1. We select only sky areas that are covered by survey 
with completeness of over 85%, so that the galaxy 
field nearly uniformly samples the underlying den- 
sity field. 

2. We compute the overdensity field and embed the 
survey volume in a periodic box (of size 1762"^ 
(ft.~i Mpc)3, with cells of size 6.9^ [h-^ Mpc)^). 
We assume that the region without observations is 
at the mean density of the universe. 

3. We compute the Wiener filtered density field and 
also the velocity field (via our six-step reconstruc- 
tion method as described in Sec. |T]). 

4. We compute the momentum field as it is projected 
along line of sight and calculate the 2D momentum 
field that would be used to cross-correlate with the 
appropriate CMB field to get the Cf^^. 

We made the templates available in both ACT and var- 
ious PLANCK resolutions. An update of the templates 
with full SDSS data-set will be made as data and time 
allow. The templates are made available to public on our 
website Their respective (^QQ-^^'yt^^ fQj. ^j^g tem- 

plates are also available as the filtering we use in the pa- 
per affect the overall normalization of the signal. In gen- 
eral, it is simpler to compare the provided (^QS>,obi,oivcd 
and the observed Cf^^ when one cross-correlates the tem- 
plates with the CMB sky. 

The two bias and gas-dependent parameters, ^off (the 
weighted visibility function) , and 6off (the weighted bias) 
are complicated, but slowly varying, functions of the 
galaxy and gas fields; we take our fiducial model {R of 
unity) to have constant, scale- independent bias b of 1.2 
(2) for main galaxy (LRG) sample, and fif, of 0.0441, h 
of 0.71 with Xe of 1.0. For reference, we plot the 2D kSZ 
field made from SDSS DR4 Main galaxy template here 
in FigEl 
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FIG. 4: Plotted is the 2D kSZ field made from SDSS DR4 
Main galaxy momentum template. The template predicts 
both hot and cold regions (due to the fact that this is a mo- 
mentum field). This implies that the template is likely to be 
nearly orthogonal to effects that are scale with the galaxy 
or the matter density (such as integrated Sachs-Wolfe ef- 
fect ([131), Rees-Sciama effect [ssl ]. tSZ, lensing, dusty galaxy 
emission and to the galaxy foregrounds. 




FIG. 5; The 2D kSZ template reconstructed from SDSS DR4 
Value-added Galaxy Catalog (VAGC.) 



V. DISCUSSION 

A. Possible Caveats 

Our method deals well with a number of possible in- 
terfering signals; we discuss them here. 

1. Thermal Sunyaev-Zel'dovich. The tSZ contributes 
to our estimates as a noise term, since tSZ sig- 
nals are uncorrelated with the velocity fields; our 
momentum-kSZ correlation is not biased by the 
presence of tSZ. In principle, we can also use mul- 
tifrequency data from CMB to separate the tSZ 
signal from kSZ as tSZ and kSZ have different spec- 
tral signatures. We can then reduce the noise con- 
tributed from tSZ. There is in general leakage of 
the tSZ signal into CMB maps due to uncertain- 
ties in the detector frequency response function and 
due to relativistic tSZ effects (see Ref. ^^). Other 
contaminants that can not be picked out with fre- 
quency information, such as the Rees-Sciama (non- 
linear ISW) and lensing effects, are generally lower 



in amplitude; they, too, contribute to the noise, but 
again in an unbiased fashion. 

2. Bias. Bias is treated as redshift and scale depen- 
dent throughout our theoretical calculations. How- 
ever, as the ratio between the observed and theo- 
retical Ci are degenerate between the effective bias 
(6^ff ) and f/off , one can determine bias (with its red- 
shift and scale dependence) but not being able to 
learn very much about the missing baryons. 

Since one can calculate bias from the power spec- 
trum analysis of galaxy populations, we assume 
that bias can be determined fully in the large scale 
structure survey before one uses the survey to pro- 
duce the gas-momentum field of the Universe. 

For simplicity in our template production we as- 
sume scale and redshift independence for bias 
(which is appropriate for the scale and redshift we 
are working with) in the production of the tem- 
plates and take on values of bias which are prede- 
termined via other efforts. 

In particular we use bias determined by [2TI| in the 
construction of the filtering functions Wgai and Wy . 
In the calculation of signal-to-noise, we assume a 
bias for the galaxy population we (will) use to con- 
struct the template for various experiments. 

3. Dusty Galaxies. Star forming galaxies can be very 
luminous at submillimcter wavelength range, as 
suggested by Submillimcter Common User Bolome- 
ter Array (SCUBA) and other experiments (see [s^ 
for a review). The emission is mainly from dust 
around the galaxies at approximately '~ 10 K, ra- 
diating mainly at sub-mm wavelength, but can also 
be significant at millimeter wavelengths. Fortu- 
nately, the dusty galaxies do not correlate with the 
bulk flows of the Universe and thus, by the same 
argument presented for tSZ, do not bias the mea- 
surement of the gas-momentum-kSZ correlation. 

4. Galactic Foregrounds. Galactic Foregrounds affect 
mostly the large-scale of the CMB and the galaxy 
observations and since they do not correlate with 
the large-scale velocity field, we can safely assume 
that it would not bias our estimator, but would at 
most be adding to the noise term. 

We could also consider reducing the contribution of 
the noise from the galactic foregrounds in our cross 
correlation by limiting ourselves to certain multi- 
poles where the galactic foregrounds are weak. We 
can gauge the effect of the galactic foregrounds by 
cross-correlating the momentum field with (for ex- 
ample) dust extinction maps (60| to find out the 
range of multipoles that are not affected by the 
galactic extinction. 
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B. Further applications 

We can use the gas-momentum-kSZ temperature cor- 
relation to constrain not only the baryon content of se- 
lected regions of the Universe, but also the baryon profiles 
and evolution of its density. We give examples below to 
describe how in theory we can apply the gas-momentum- 
kSZ temperature correlation to understand gas profiles 
and baryon evolution over redshift: 

1. Gas profiles: Assuming that we know the bias and 
ionization fraction of the Universe in the same re- 
gion of interest, we can parameterize the electron 
density as a function of radius (from the center of 
the nearest galaxy or region under consideration) 
and generate the appropriate visibility functions. 
We can then construct momentum templates and 
maximize the cross-correlation coefficient between 
the electron density and the CMB to determine the 
best-fit parameters. This can be applied to finding 
gas profiles in voids or around different types of 
galaxies. 

2. Baryon content at different redshifts. With ever 
increasing size and volume coverage of galaxy red- 
shift surveys, we can create momentum templates 
based on galaxy surveys binned at different red- 
shifts. We can then trace the evolution of electron 
density through a large redshift range. 

VI. CONCLUSION 

We have presented a new method of generating a tem- 
plate for the kinetic Sunyaev-Zel'dovich that can be used 
to detect the missing baryons. This momentum template 
is based on computing the product of a reconstructed ve- 
locity field with a galaxy field. Since the kinetic Sunyaev- 
Zel'dovich effect is a line of sight integral of electron mo- 
menta (modulo constants such as Thompson scattering 
cross-section), the combination of a galaxy redshift sur- 
vey and a CMB survey can detect the ionized gas. 

Unlike other techniques that look for hot gas or met- 
als, this approach directly detects the electrons in the 
IGM through their signature on the CMB. Since the kSZ 
is produced simply by Thompson scattering from free 
electrons, there is no need for detailed knowledge about 
the metallicity of the medium and its evolution, nor for 
an understanding of potentially out-of-equilibrium level 
populations. 

The kSZ amplitude is the product of the gas density 
and its velocity, and the latter, as we have shown, can be 
well-modeled by linear theory, providing an independent 
check, on different scales and relying on different physics, 
from baryon surveys that employ tracers that scale as the 
square of density, or that are more concentrated in the 
center of the most massive virialized objects. 

Our studies in this paper find that a CMB survey, with 
sufficient resolution to push past the "damping tail" of 



the primordial fluctuations but not necessarily covering 
the entire sky, partnered with an overlapping galaxy sur- 
vey with sufficient number density to find structures at 
these angular scales, can provide an unambiguous detec- 
tion of the baryons. 

We have estimate the expected signal-to-noise for de- 
tecting the galaxy-momentum kSZ cross-correlation for a 
few different combinations of current and upcoming ex- 
periments. The estimated signal-to-noise for detecting 
the galaxy-momentum kSZ cross-correlation is 4.1, 8.5, 
12.2 for ACT (with survey area of 2000 deg^) with SDSS- 
DR4, SDSS3 and ADEPT. The estimated signal-to-noise 
for PLANCK with SDSS-DR4, SDSS3 and ADEPT is 
11.2, 22.5 and 32.2. The proposed estimator provides an 
exciting avenue into understanding the ionized gas in the 
Universe in the near future. 

We have produced momentum templates from Sloan 
Digital Sky Survey, so that the scientific community make 
take arcminute scale CMB data from current and upcom- 
ing experiments such as ACT, SPT and PLANCK and 
determine not only the gas fraction of the Universe, but 
also the distribution and evolution of the gas fractions 
and gas profiles around different regions in the sky. 
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